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Abstract 

We analyze the problem of electronic transmission through semi-infinite and finite regions of graphene which are characterized 
by different types of connections between the Dirac points. These valley symmetry breaking Hamiltonians might arise from 
electronic self-interaction mediated by the dielectric environment of distinct parts of the substrate on which the graphene sheet 
is placed. We show that it is possible to have situations in which we can use these regions to select or filter states of one desired 
chirality. 


I. INTRODUCTION 

Graphene is a truly two-dimensional (2D) crystal with 
Dirac-like quasi-particles mm- These charge carriers 
are the result of the wavefunction interference due to the 
honeycomb lattice structure and hence they cannot exist 
outside the many-body system (electrons-plus-lattice). 
From this perspective, the low energy, long wavelength, 
Lorentz invariance in this system is an emergent phe¬ 
nomenon and, as such, can be influenced by external 
factors such as disorder, applied electric and magnetic 
fields, and structural deformations of the lattice such as 
pressure, strain, and shear |3]. 

The low energy physics of these Dirac quasi-particles is 
strongly influenced by the fact that the two Dirac cones 
sit at the corners of the hexagonal Brillouin zone (K and 
K’). These cones are related by time reversal symmetry 
and define the chirality of these quasi-particles in terms of 
their momentum relative to each cone. Unlike the case 
of neutrinos, where only one chiral flavor exists, Dirac 
particles in graphene have both flavors, and they can be 
either right-handed or left-handed depending on whether 
they reside in one cone or the other. While in a per¬ 
fect, non-interacting, graphene sheet these chiral states 
are decoupled, in a real graphene sample they can be cou¬ 
pled by the external perturbations mentioned above. At 
low energies, one can have either intra-valley processes, 
with small momentum transfer, that preserve chirality, 
or inter-valley processes, with large momentum transfer, 
that mixes chiralities. It is well-known that, in the pres¬ 
ence of weak disorder, intra-valley processes lead to weak 
anti-localization effects |3] , whereas inter-valley processes 
lead to weak localization (5]. Hence, the coupling between 
K and K’ points play an important role in the physics of 
graphene. 

In this work we study the transport through regions 
where the K and K’ points are not coupled and regions in 
which they are actually coupled. Although the Hamilto¬ 
nians with which we deal are manifestly non interacting, 
these couplings might arise from the interactions between 
the electrons, as we argue below. Our main goal is to un¬ 
derstand how the Dirac electrons behave across the in¬ 
terfaces separating those two regions as a way to classify 
the possible scattering mechanisms in graphene. Such 


a situation can also be artificially created by depositing 
graphene across substrates with different dielectric con¬ 
stants. For instance, there are quantum Monte Carlo [7] 
calculations that indicate that suspended graphene (i.e., 
on vacuum with dielectric constant Eq = 1) is a excitonic 
insulator and graphene on Si02 (esi 02 ~ 3) should be a 
semi-metal [5] . An interface between those regions would 
have a transistor-like effect with a large on-off ratio for 
current flow. 

Our starting point is the well-known low energy effec¬ 
tive Hamiltonian for neutral graphene that is given by [T] 
(we use units such that vp = h = 1): 



where 'F| = ^aj,6|^, aj and bj being the creation op¬ 
erators for electrons in sublattice A or B, respectively, 
in the Dirac cone i {i = 1,2) and a = {ax,ay) is a 2D 
vector whose components are Pauli matrices (we ignore 
spin variables). 

In the representation of the Hamiltonian Q there are 
many ways to couple the K and K’ points [5]. These 
couplings have different symmetries and hence represent 
different physical processes. Notice, however, that all the 
processes that couple the two cones can be represented in 
terms of combinations of the identity matrix and Pauli 
matrices in the "valley" space. We will deal only with 
real coupling potentials in such a way that we can write 
explicitly all the generators of the different couplings as 
combinations of basic matrices: 

and 

v,=a(';j), (5) 
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where A is a parameter that is assumed to be positive 
throughout the paper since negative deltas will bring 
no new physics. Of course the most arbitrary poten¬ 
tial would be a linear combination of these with different 
delta parameters (and in this case the sign difference be¬ 
tween them might be relevant). 

As we are going to show, only Vy and Vx are of real 
interest and will be the focus of our studies. We will also 
show that the transmission through regions described by 
these matrices presents very unusual properties and can 
presumably be measured experimentally (we stress, nev¬ 
ertheless, that it is not our aim to fully solve the transport 
problem considering finite size and disorder effects, but 
to analyze the physics of our proposed model which we 
believe can account for the main properties behind, for 
instance, the contact resistance of nanoscopic graphene 
junctions). 

The paper is organized as follows. In Sec. II we give a 
very brief introduction to the problem of "free" graphene 
and some of our notation. In Sec. Ill we present and 
analyze the Hamiltonians with which we work, giving a 
somewhat general notion on why our choices of interac¬ 
tions and their nature, solving the corresponding eigen¬ 
value problems and analyzing the symmetries they obey. 
In Sec. IV, we (less) briefly discuss the transmission prob¬ 
lem in two dimensions. In Sec. V we solve the problem 
of transmission through the interface between two semi¬ 
infinite regions of graphene, one with "free" electrons and 
the other with self-interacting quasi-particles and finally, 
in Sec. VI, we solve the finite-barrier-like problem of 
transmission through a finite region of the interacting 
material. In Sec. VII, we try to make a connection of 
our results and the Landauer formalism. In Sec. VIII we 
present our conclusions. 


II. BASIC PROPERTIES 


The eigenvalue problem given by 0 can be solved by 
going into the momentum representation 


with energy eigenvalues, E = ±fc, and eigenstates of well 
defined momenta around the Dirac points (with their cor¬ 
responding spinors, see below). Going back to the posi¬ 
tion representation one has V'K(K')(r) = V’±K(K')(k)e*'^ '’ 
where the spinors are 


V^±K(k) 


V'±K'(k) 


1 


1 
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with 9 = aicta,TL[ky/kx) and the vector k, whose mod¬ 
ulus is connected to the Fermi energy by the disper¬ 
sion relation, is centered either around K or K', respec¬ 
tively. These eigenstates, which are symmetric and anti¬ 
symmetric linear combinations of the states referring to 
the A and B sublattices, allow us to introduce a particle- 
hole representation which reads 


^ A: fljk,i - 


2 = 1,2 


( 8 ) 


where the hole states have negative energy. Thus, we de¬ 
scribe the electron behavior in the free graphene by a the¬ 
ory of non-interacting massless fermionic quasi-particles 
with two different "flavors" (and their corresponding 
anti-particles). 


III. INTERACTING HAMILTONIANS 

Electron-electron interaction plays an important role 
in graphene physics and transport problems [inilll]. Its 
effects depend strongly on the dielectric function and is 
independent of the electronic density m- Since the di¬ 
electric function can be affected in many ways, as, for 
instance, by adsorbed atoms on the surface or deposi¬ 
tion substrates, we propose that there may be situations 
in which the resulting effective electron-electron interac¬ 
tion can cause mixing of the Dirac cones and we show 
how it could happen, at least in a heuristic mean-field 
approximation. Our approach is phenomenological, and, 
therefore, the specific shape of the interaction potentials 
will not be given and should be chosen from microscopic 
arguments. There are many possible situations, and one 
of them is the possibility that the (Coulomb) interaction 
becomes screened [innn. Hence, we justify the choice of 
our approach and will make the symmetry assumptions 
we find reasonable to achieve our results. 

So, let us suppose that the electron-electron interaction 
contribution to the usual tight binding Hamiltonian is of 
the form 


Hj — ^ ^ UijTliTlj . (9) 

Separating the terms in the summation corresponding to 
each sublattice, we arrive at 

Hia = ^ UijUAinAj ( 10 ) 

Hjb = ^ UijUBinBj ( 11 ) 

Hjab — ^ ^ hijBAi^Bj ■ 
i,3 


( 7 ) 


2 


We have assumed that the interaction between the elec¬ 
trons within each sublattice is equal, which actually need 
not be. 








Now we perform the mean field approximation. The 
values of the mean-fields should actually be determined 
in a self-consistent way, which will not be done here, and 
be justified by microscopic arguments. Our point is that 
if there is any physical process which gives rise to these 
mean-fields, the proposed decoupling will lead to our re¬ 
sults. We make more comments on the physical nature 
of this approximation further below. That being said, we 
can choose, for example, the first Hamiltonian and get, 
to first order in the deviations from a mean density of 
electrons in sublattice A, 

Hi A « ^ [uAi (riAj) + nAj {nA%) - (nAi) (nAj)] 

= ^iUAi — const., (13) 

i 

where A is the effective on-site potential generated by the 
self-interaction between the particles. In Fourier space, 
we can write 

(«) 

where Ns is the number of sites on the lattice. Hence, 
we have, 

HjA ~ ; ( 1 ^) 

kk' 

where we have ignored the constant term. Therefore, this 
Hamiltonian couples the k vectors . If we keep only the 
the long wavelength contribution to the expression just 
above, we find that this term leads to the four couplings 
(following our notation, ok —>■ ai,aK' —>■ 02 ) 

a\ai; a\a 2 (16) 

and 

a\a2; a|ai. (17) 

The first two lead to diagonal terms which are of no in¬ 
terest to us. The two other give contributions exactly 
of the type we want. The same arguments can be car¬ 
ried over to Hjb- Since the potential must be such that 
Uij = Uji, its Fourier transform is real. The mean fields 
should also be real and this justifies our choice of only 
analyzing the four real V Hamiltonians. As for the phys¬ 
ical realization of these mean fields, we believe that the 
presence of the substrate might, in specific cases, affect 
the electrons in the material and break the symmetries 
between the sub-lattices in such a way that ordered states 
like, for example, a CDW could develop. More thought 
in this direction is needed but will not be pursued in this 
paper. 

If one wants a more thorough approach, one can use 
the Alicea and Fisher model m- There it is shown how 
the long and short range interactions give rise to terms 
like staggered densities in each sublattice. 



Figure 1: (color on-line) Schematics of the dispersion relation 
for the Hamiltonian in the presence of the Vy term (arbitrary 
units). 


where UA.Bi, i = 1 , 2 , refers to the number operator in 
the sublattice A or B and Dirac cone 1 or 2. In the Alicea 
and Fisher model, the long wavelength Hamiltonian has 
(local) contributions due to electronic densities given by 

Ptot + pltag = a\aia\ai + 4020^02 -b 

+b\bib{h + blb2b\b2 + b\biblb2 + 
+a\a 2 a\ai + a\aiala 2 + 62 ^ 24 ^ 1 - 

(19) 

This term comes from the on-site repulsion interaction 
(a more detailed analisys can be found in the original pa¬ 
per m)- We can see that a mean-field decoupling might 
also generate our Hamiltonians in this case, however, the 
mean-fields in this case result directly from the replace¬ 
ment 0102 —f ^Oia|^ , for example. We find this type of 
mean-field harder to justify physically but leave it here 
as an example. 

As we had formerly proposed, we have given some 
heuristic arguments for our ad hoc choice of the form of 
the Vi interactions. We proceed now to exploit the con¬ 
sequences of these model Hamiltonians in some specific 
cases. 


Spectrum in the presence of Vy 


In the presence of the disturbance ^ 
is given by 


the Hamiltonian 
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( 20 ) 


The solution of the eigenvalue problem gives Ea^b± = 
±\/fc2 -b = ±e (Fig. 1) with a gap of size A in the 
spectrum. The spinorial part of the eigenvectors (we omit 


Pstag = tlAl + nA2 — n_Bl — ^32, 


(18) 
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the plane waves for simplicity in this whole section) 
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( 22 ) 

where A± = ^ Here, A and B do not refer to the 
different triangular sub-lattices but to the different de¬ 
generate states originating from the dispersion relation. 

as we saw, can use eq.Q to put this Hamiltonian in 
a particle-hole representation, where particle and hole 
states are described as symmetric and anti-symmetric 
combinations of the different sub-lattices wave functions. 
The effective interaction between the particles and holes 
is then written as 


^ph — 


/ k 0 0 

0 -k -A 
0 -A fc 
\ A 0 0 



(23) 


where ph refers to the particle-hole representation. Eq. 
(201 can be seen to obey orthogonal time reversal sym¬ 
metry (exchanging the valleys by time reversal) and it 
is seen to generate a band gap. We see that this kind 
of effective potential is equivalent to the introduction of 
an asymmetric coupling between particles and holes from 
different Dirac cones. 


Spectrum in the presence of Vx 


In the presence of a perturbation of the form (|^, the 
full Hamiltonian becomes : 
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(24) 


Diagonalization will lead us to E±± = ±k ± A (Fig. 
2) with no gap in the spectrum but a shift of the Dirac 
cones relative to each other by A. The eigenvectors are 
now 
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Figure 2: (color on-line) Schematics of the dispersion relation 
for the Hamiltonian in the presence of the 14 term(arbitrary 
units). 
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, (26) 


and, in the particle-hole representation we have 


Hph — 


fc 0 A 0 \ 
0 -fc 0 -A I 
A 0 fc 0 

0 -A 0 -k J 


(27) 


The Hamiltonian now obeys S time reversal symmetry 
and there is no energy gap (5 is the sympletic time 
reversal symmetry in the sense described in m)- 
However, there is another interesting property that 
shows up. Because of the different translations in 
energy of the two Dirac cones, there can be phenomena 
akin to the Klein paradox in transmission problems. 
The particle-hole representation reveals an effective 
interaction that favors the coupling of different valley 
particles over the coupling of different valley holes states. 


Spectrum in the presence of Vi and Vz 


Two other Hamiltonians that possibly connect valleys 
are given by Q and (|^ and therefore. 
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and 
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These are of no actual interest to transmission prob¬ 
lems. The dispersion relations will be, respectively, 

E±±^ = ±A/fc2-f A2±2fcAcos0 (30) 


and 


E±±y = ±\/fc2-hA2±2fcAsin6». (31) 

These problems are related to each other by rotations of 
7r/2 around the z axis, and hence do not introduce any 
new physics. Besides, we see that the effect of this po¬ 
tential on the dispersion relations of the particles reduces 


to just a change by ±A in the x (first case) or y (second 
case) components of the particles’ momenta, which can 
be renormalized and the resulting quasi-particles behave 
no differently from free particles. Therefore, in the re¬ 
gions affected by these potentials, the transmission prob¬ 
abilities will be equal to one, as can be shown by solving 
the semi-infinite or finite interacting regions problem. We 
will not deal with these two cases any further and show, 
only for completeness, the expressions for the spinors as¬ 
sociated with the Hamiltonians above, which respectively 
read: 
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IV. THE TRANSMISSION PROBLEM IN 2D. 
"BARRIERS" AND "STEPS" 

From now onwards, we will be dealing with situations 
in which there is transmission of electronic waves from 
"free" electron graphene into semi-infinite or finite re¬ 
gions of "disturbed" graphene, and also reflection back 
into the original region. We call these two situations step 
and barrier problems, respectively. 

Although there has been extensive use of wavefunction 
matching to describe problems of electronic behavior in 
graphene (some examples can be seen in |18| and |19|h we 
would like to make some comments of our own. We shall 


start by briefly talking about the simple 2D scattering 
from straight interfaces. Although it is a simple problem 
it is very useful to establish the terminology we employ 
in the more complex cases. All problems of transmission 
begin with the determination of the wave functions in 
different regions and the enforcement of the boundary 
conditions they obey. We assume that the particle in 
medium I has positive energy, with momentum on the 
Dirac cone K and moves to the right. Since both types 
of interactions connect states of the different Dirac cones, 
conservation of the valley "flavor" need not take place any 
more and the states corresponding to medium I must be 
given by 
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(34) 


Moreover, if the problem is of a barrier type, we should particles in medium III must also have positive energy 
be able to find in medium III the particles with mo¬ 
menta around both valleys, independently of its flavor 
in medium I. Conservation of energy demands that these 
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and then we have 
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The states accessible to the particles in medium II depend 
on the dispersion relations of the Hamiltonian and we will 
deal with each specific case in the next sections. 

Since the Hamiltonians are of first order in P, we need 
only to invoke the continuity of the wavefunctions at the 
interfaces of the different media to be able to determine 
the coefficients and However, some care must be 
taken when evaluating the transmission and reflection 
amplitudes in the step problems. These must be deter¬ 
mined by the conservation of the probability current in 
the direction normal to the interface (in our case, the x 
direction). Beginning with the continuity equation we 
have, for medium i = I, H, HI, 

= (36) 

which in the stationary regime implies that pi = |'I'i|^ is 
time independent and we have 



Region I Region II 


Figure 3: (color on-line) Schematics of the step transmission 
problem for Vy Hamiltonian. The red (light gray) curves are 
associated with hole states and the blue (dark gray) ones with 
particles states. 


must be used to find the correct expression. The x com¬ 
ponent of the probability currents in media I and HI will 
be always the same in our problems and are found to be 


Jlx 


l-kil" 



and 


JlIIx = 




cos 9. 


(41) 


(42) 


V • Ji = 0. (37) 

Now, since our problem is translational invariant along 
the y direction, is independent of the y variable and it 
finally reduces to a one dimensional conservation problem 

Jix = const. (38) 

Although the conservation of the probability current is 
relevant only in one dimension, the coefficients Vi and 
ti and the eigenstates (with which we will calculate the 
probability current) all depend on the angles that the 
momenta of the particles make with the normal to the 
interfaces in each medium. Since the states are equal to 
each other at the interfaces (boundary conditions), we 
have that 

Jlx — Jllx — JlIIx: (^^) 

the last equality happening only for barrier problems. 

We calculate the probability current in each medium 
by taking the mean values of the current operator defined 

by 

£TJ 

J = -- = e.^, (40) 

with the states 4*^, where A is the electromagnetic vector 
potential. Notice that depending on the type of interac¬ 
tion (for instance, if it is momentum dependent), the ex¬ 
pression for the current will not assume this usual simple 
form and the definition through the functional derivative 


With these, we see that, for barrier problems 

Jlx = Jiiix => 1 = |ci|^ -b |r2|^ -b |ti|^ -b 1 ^ 2 !^ , (43) 

and we recover the usual result. In step problems, we 
will have different coefficients, depending on the states 
involved in the scattering at the interface. 


V. SEMI-INFINITE REGIONS: STEP PROBLEM 

Vy interaction 

This problem is very much alike the usual quantum 
mechanical problem of transmission through a potential 
energy step smaller than the particle’s energy. In Fig. 
3, we show the energetics of the problem and we easily 
see that for k < A, k being the energy of the incident 
particle, there will be no transmission, and hence we need 
the energy of the incident particle to be larger than A. 
Moreover, conservation of energy and momentum along 
the y direction demands that 

k = \/~q^~+~A^ (44) 

ksinO = qsina. (45) 

where a is the angle between the y and x components of 
the wavevector of the transmitted wave. The transmis¬ 
sion can take place with any of the two degenerate states 
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of positive energy, and, consequently, the boundary con¬ 
ditions at the interface leads us to the system 
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(46) 


Its solutions are 


ri = 

q sin a — k sin 9 

(47) 

k + qcos{a + 9) ’ 

r2 = 

. cos 9 

lA 

k + qcos(a + 9) ’ 

(48) 

tl = 

,—^-- cos 9 cos (^^) 

72 fc fc + g , ; 2 J . 

' k + qcos{a + 9) 

(49) 

^2 = 

i\/2k{k - q)- ) ^ 
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which shows that there are reflection and transmission to 
all the energy degenerate states. Similar expressions for 
the transmission and reflection amplitudes and probabil¬ 
ities (which follow below) are always found when dealing 
with wave-function matching in graphene and have pre¬ 
viously been derived in the literature UHlIISl. 

The evaluation of the probability current in medium 

II reveals that Jjix = f + 1^21^^ cos a. Hence, by 
([ 3 ^ and Cil, 


T=(\t 




q COS a 
k cosO 


(51) 


We can use (^44|, the coefficients (47l-(50) and the defi¬ 


nition of a = {qy/Qx) , to express this transmission 

probability (TP) as a function of the given parameters of 
the disturbing potential and of the incident wave. 

So we get 


T = 


2 -\/l — coso; 

1 -I- y/l — 77 ^ cos a ’ 


(52) 


with the angle of emergence given by 


a = tan 


sin 6* 


-y/cos^^^^^ 


(53) 


where we defined r] = A/k (this procedure will be re¬ 
peated throughout the whole paper in a more succinct 
fashion, unless some special warning is necessary). 

One should be careful about some details of these so¬ 


lutions. First of all, one notices from (531 that there 


are situations in which, even if the energy of the particle 
is greater than the interaction energy, for large enough 
incidence angles, a will be complex, also leading to a 
complex TP, which has no physical meaning. This phe¬ 
nomenon is analogous to total reflection in electromag¬ 
netism. This happens because although in the transmis¬ 
sion and reflection probabilities (TP and RP) expressions 
the transmission and reflection amplitudes appear only in 
squared moduli, the probability current also depends on 
the angles. Making the substitution of imaginary qx and 
a first in the wavefunctions and then proceeding with 
the calculation of Jux reveals that actually Jux = 0, 
as expected. This is the same as happens for scattering 
by step potentials in usual quantum mechanics when we 
deal with energies above and below the step energy. 

In Fig. 4, we show the general behavior of the TP. 
Some interesting features of the problem arise now, as we 
see that state 2 only contributes to transmission about 
the extreme values of the allowed angles. As we raise the 
barrier in relation to the particle’s energy, we see that the 
probability of transmission drops to zero, as expected. 
We also note that there is a localization of the beam, 
with respect to the possible incidence angles which allow 
for transmission with non-zero probability. Looking at 
the behavior of the "angles of refraction", which charac¬ 
terize the direction of propagation of the wave-fronts of 
the spatial parts of the wavefunction, as a function of the 
angle of incidence, we see that for small values they are 
about the same, at least for low rj. One should notice 
that these angles need not have actual relation to the 
direction of propagation of the probability current, since 
the spinors also depend on the momenta and this affects 
the mean value of the current operator (it is connected 
to the spinors through the a matrices). 


Vx interaction 


The scheme of the dispersion relations for this physical 
situation is shown in Fig. 5. 

We notice that there are two non-degenerate bands ac¬ 
cessible to the incident particles (we will call these the 
states 1 and 2), which lead to particles moving with dif¬ 
ferent momenta inside the medium, due to the conserva¬ 
tion of energy. This will lead to two different angles 
i = 1, 2, for the emergent particles on the right side of 
the interface, related to the angle of incidence 9 by 

g^sinOi = fcsin0. (54) 


There can also be two different cases, which will change 
the states accessible to the emergent particles, namely 
k > A and k < A. These considerations lead to the 


following situations: 

i) k > A In this case (Fig. 5 left), the accessible 

states in medium II are given by |x-i-+) and |y_|_) from 

equation (25 \ and we have 


k = qi+A = q2 — A (55) 
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Figure 4: (color on-line) Behavior of the TP for the Vy step problem. In the upper and lower figures on the right, we have 
rj equals to 0.3, 0.6 and 0.9,for the green (light gray), red (mid gray) and blue (dark gray) curves, respectively. (Upper left) 
Total TP as a function of the angle of incidence. (Upper right) We separate here the TP associated with states 1 (full line) 
and 2(dashed line). (Lower left) TP for normal Incidence as a function of the ratio of the energies of interaction and incidence. 
It’s important to notice that the contributions here come only from state 1. (Lower right) "Refraction angle" as a function of 
the incident angle for the same values of 77 above. 



Region I Region II Region I Region II 


Figure 5; (color on-line) Scheme of the dispersion relation for the step transmission problem for Hamiltonian Vx- The red (light 
gray) curves are associated with hole states and the blue (dark gray) ones with particles states, (left) Transmission for k > A. 
(right) Transmission for k < A. 


and the boundary condition at the interface leads to 
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The solution of the linear system in this case is 


The probability current in the x direction is given by 
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r2 
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1 [sin (_ cos ) sin 9] 
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(58) 
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Jiix = e cosai -I- 1 ^ 2 !^ cos 02 ) , 


leading to the transmission coefficient 
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Substituting the transmission amplitudes of the solution 
gives 


T = cos 6 


cosai 


cos a 2 


1 + cos[ai -\- 0] 1 + cos[a 2 + 0] 


(63) 


Proceeding in the same way as in the last section, we get 


ai 


02 


tan ^ 


tan ^ 



(64) 


(65) 


II "less refractive" and, consequently, are fully reflected. 
However, for particles of band 2, the medium is "more 
refractive" and there will always be transmission. 

We notice that the main contribution to the total TP of 
states of band 1 (those subject to total reflection) comes 
mostly from small angles 6, so the current must be al¬ 
most normal to the interface if one expects to transmit 
these particles. Looking at the "angles of refraction" as 
a function of the angle of incidence, one can easily rec¬ 
ognize which band is subject to total reflection. 

ii) k < A In this case, we have the possibility of 
conversion of a particle state into a hole state, as shown 
in the right diagram of Fig. 5, and the conservation of 
energy gives us 


with which, one can plot the behavior of the TP, analo¬ 
gously to what we have shown in Fig. 6. 

Conservation of the y component of the particle’s mo¬ 
mentum would reveal that, in analogy to the classical 
electromagnetic case , particles of band 1 feel the medium 


k = —qi + A = q 2 — A. 


( 66 ) 


The linear system to be solved involves now the states 

lx_h) and Ixh _), from (251 and (261 respectively, and 

reads 
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Its solution is 


ri= - 


1 [cos (■ 


sin 6] 


sin cos 


r2 ■■ 


1 cos(^ii^)cos6» 
2 sin(^i!^) cos 


a2+6 ^ 


( 68 ) 


(69) 


U = i 


cos 9 




(70) 


t2 = 


COS 6 


\/2 sin ( 


Q 2 +^ ^ ’ 


and the probability current and TP are given by 

Jiix = e |<i|^ cosoi-I- 1 ^ 2 !^ cosa 2 ) , 
cosai |, ,2 , cos 02 I ,2 


cos 9 


(71) 

(72) 

(73) 


which, after proper substitution of the coefflcients leads 
to 


T = cos 9 < — 


cosai 


cos 02 


1 — cos[ai 9] 1 -|- cos[a 2 -I- 9] 


■ (74) 


One can see from this expression the expected appear¬ 
ance of the Klein paradox, by looking at the denominator 
of this expression’s first term in the case of normal inci¬ 
dence. The solution of this is the well-known argument 
that the hole-like particle must have an inverted momen¬ 
tum, so that it will continue traveling from left to right 
in the expected direction [5]. For such, we demand 



We notice that, with these conventions for the angles, the 
expression (631 is actually valid independently of the en¬ 
ergy of the particles and their interaction energy, but we 
will keep the solutions separated. The general behavior 
of the TP for this case is shown in Fig. 7. 

We notice that as we keep raising the barrier (or lower¬ 
ing the particle’s energy), the contributions of both bands 
tend to become equal. It is interesting to see that when 
rj = 2, particles of band 1 go through the barrier as if 
there was nothing there, independently of the angle of 
incidence. The behavior of the particle in band 1 as a 
function of 77 > 1 is such that the transmission starts 
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Figure 6: (color on-line) Behavior of the TP for the Vx k > A step problem. The curves plotted are for r] equals to 0.3, 0.6 
and 0.9 corresponding to green (light gray), red (mid gray) and blue (dark gray) respectively. (Left) Total TP in function of 
the angle of incidence. (Middle) Contributions to the TP given by the separate states 1 (full line) and 2(dashed line). (Right) 
"Refraction angles" as a function of the incidence angle for the same values of rj above. Full lines correspond to a\ and dashed 
lines to Q 2 . 


around zero for ry « 1 increases to uniform probability 
when rj = 2 and then tends to the same type of behavior 
of the particle in band 2. 

From the behavior of the angles ai and 02 as a function 
of 9 it is clear how, as we keep raising the barrier (or 
lowering the energies), the direction of propagation of the 
wave-fronts tend to become equal. For rj = 2 there is the 
explicit inversion of behavior for the particle in band 1 
and it stops being subject to total reflection (notice that 
the rightmost plot of Fig. 7 is actually of ai — tt but, 
since the holes move contrary to their momenta, this is 
the actual direction of motion of the particle’s wavefront). 

We can use the results of these last two cases to pre¬ 
pare a quasi-particle current in the material in a definite 
quantum state or even in a known linear combination of 
states. For the normal incidence problem, we also see 
that there will be full transmission, independently of the 
relation of the energy of the particle to the strength of 
the interaction. 


VI. FINITE REGIONS:BARRIER PROBLEMS 

Now we revisit the problems of the last sections for 
the case of a finite region of disturbed graphene. The 
mathematical developments proceed almost in the same 
way as before, the only difference being that we need 
to match the wavefunctions from medium I and medium 
II at a: = 0 and from medium II and medium III at 
X = d, d being the width of the barrier, which will lead 
to a more complex linear system. Hence, we have the 
following situations. 


Vy interaction 

The matching of the wave function leads to the system 
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Figure 7: (color on-line) Behavior of the TP for the Vx k < A step problem. The curves plotted are for r] equals to 1.4, 2.0 
and 8.0 corresponding to green (light gray),red (mid gray) and blue (dark gray) respectively. (Left) Total TP In function of 
the angle of incidence, (middle) Contributions to the TP given by the separate states 1 (full lines) and 2(dashed lines), (right) 
behavior of ai — 7r(full lines) and a2 (dashed lines) as a function of 6. 


Notice that in medium II, we can now have both trans¬ 
mitted and reflected waves, because of the presence of 
two interfaces. We can solve this to And the coefflcients 
Ai and Bi in terms of U and then solve a problem anal¬ 
ogous to those of the last sections. The dependence of 
Jjix on the position variables can be easily seen to van¬ 
ish, as one should expect, and the resulting transmission 
probability can be directly obtained, 

T=\hf + \t2f. (79) 

The expression for the coefflcients are cumbersome, but 
we show the ones corresponding to the ti’s and r^’s. They 
are: 


the incident and emergent angles, for they are the same. 
We also have the situation of total reflection in this case 
although, due to tunneling, the kinky behavior of the TP 
we had in the step problems is not apparent any more. 

One should notice the oscillatory behavior generated 
by the self interference of the wave reflected inside the 
medium II (Fabry-Perot like interferences). In the middle 
plot, we realize that these oscillations are affected mainly 
by the width of the barrier. We also notice that raising 
the relative value of A over k causes again a focalization 
of the beam and gives rise to a resonant phenomenon of 
full TP. 

These resonant points can be easily analyzed. The TP 
for 9 = 0 reduces to 


ri 


ri 


sin {dqx') {—q sin a k sin 9) 
q cos (dqx) cos a cos 9 + i sin {dqx) [k — q sin a sin 0) ’ 
\Jk‘^ — q^ cos 9 sin {dqx) 

q cos {dqx) cos a cos 9 — i sin {dqx) {k — q sin a sin 0) ’ 


T\ 0 =o = 


A2 

— A^ cos^ {dq) ’ 


whose maxima are seen to be given by 


(84) 


^ e q cos a cos 9 

^ geos {dqx) cos a cos 6* — fsin {dqx) {k — g sin a sin 0) ’ 

t2 = 0; (80) 

We are not able to write the expressions for the TPs 
in terms of the ratio of the incident particle’s energy to 
the potential strength A anymore. In order to get an 
expression for the TPs, given only in terms of known 
variables, we now use 


dqn = riTT => = 



(85) 


which resembles the energy of massive relativistic parti¬ 
cles in a box. Expanding ti{9 = 0) for g around mrjd, 
we also get 


h{9 = 0)^±i-^ 


rn/2 
kn + 


( 86 ) 


a 


Qx 




(81) 


where the k'^s are, as above, the energy levels of virtual 
bound states whose inverse life-time (decay rates) are 
given by 


(82) 





= 2 1-k 



(87) 


finding 

T= , _ 

k"^ cos^ 0 — A^ cos^ {d\/k^ cos^ 9 — A^) 

The behavior of the TP is shown in Fig. 8. Since the 
emergent particles come out in a medium equal to the 
one from which they enter, we do not need to consider 


{k"^ cos^ 9 — A^) 




Similar life-times should appear for different values of the 
incident angle. We see that if, for fixed n and A, the bar¬ 
rier width gets larger, the particle life-time also increases. 
This is an expected result since the mean number of times 
the particle should reflect back and forth does not change 
and neither does its velocity, so the time spent inside the 
barrier should be longer if it is larger. If the interaction 
is stronger, i.e., the barrier is "higher", one expect the 
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particle’s quantum behavior to be more important and 
therefore, even if the particle’s energy is above the bar¬ 
rier, it will spend "more time" inside it and we expect a 
longer life-time. 

If in the normal incidence TP (841, both the interaction 
A and the energy k are small, we have, np to second order 
in and A^, 


T\e=o 


A2 


( 88 ) 


fc2- A2-hA2d2 (/c2- A2)' 
Decreasing the energy to arbitrarily low valnes leads to 


checked), it is possible to have transmission for arbitrar¬ 
ily low energies, even if smaller than A. This is a phe¬ 
nomenon analogons to tnnneling in nsual quantum me¬ 
chanics barrier problems. The wavefnnction component 
dependent on the x direction becomes a real exponen¬ 
tial, which can in some sitnations penetrate in the other 
medium. 

Notice that the transmission amplitnde for states 2 (^ 2 ) 
is identically zero. This means that there can be no flip¬ 
ping of the psendo-spin dne to the Vy coupling Hamilto¬ 
nian. The situation is different for the Vx interaction, as 
will be shown in the next sections. 


T\e^o 


1 

1 -k d2A2 ■ 


(89) Vx interaction 


Hence, if the interaction energy is low enough (or even i)k > A The linear system in this particnlar case is 
if the barrier width is small enough, as can be easily now 
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The expressions for the solutions in this case are even 
more cnmbersome than in the previons case, although 
they still resemble the expressions presented in the last 
section for this particnlar sitnation. We will show only 
the different behavior of the TP patterns related to them 
(Fig. 9). The angles ai and 0:2 and the x components 
of the momenta of states 1 and 2 are the same of those 
of the corresponding case of the previons section. We see 
that, for small angles, the resonant phenomena induced 
by the finite barrier is irrelevant. There is some selection 


of states for larger angles and only states of type 2 con¬ 
tribute to the central region. For the normal incidence, 
we also have a perfectly resonant behavior between the 
two TPs (which canse the total transmission be always 
equal to 1), as a function of the width of the barrier, 
which can be nsed to select the different states. There 
is no resonance for fixed width and variable energies of 
incident particles. 
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Figure 8: (color on-line) Behavior of the TP for the Vy barrier problem. Since t 2 = 0, we only plot the total TP. On the left 
we have d = 2.0 and A = 3.0 fixed and vary the energy k of the particles as 3.35 (gray, dashed light gray), 4.3 (green, full line 
light gray), 5.0 (red, mid gray) and 10.0 (blue, dark gray). In the middle, we fix A = 4.45 and k = 9.0 and vary the width d 
through 1.0(dashed green, light gray) , 2.0(dashed red, mid gray) and 3.0(full line blue, dark gray). On the right, we plot the 
normal TP as a function of k for d = 1.5 and A equals to 1.0 (green, light gray), 2.0 (red, mid gray), 3.0 (blue, dark gray). 



ii Xy 

11 * V # / 

Vv if* li 
W 0 1 II 

IW \\i I 

\ 1/i i/I 

11 1 v / 
Ilia/ 

\ W I 11 

\ B 1 1 1 

1 1 / IFm 

' f i 

Hiv/ Vv 




k d 


Figure 9: (color online) Behavior of the TP for the Vx k > A barrier problem, (upper left) We have d = 3.0 and A = 1.0 fixed 
and vary the energy k of the particles as 2.0 (green, light gray), 4.0 (red, mid gray) and 9.0 (blue, dark gray), (upper right) We 
separate here the contributions to the total TP from the states l(full lines) and 2 (dashed lines). (Lower left) Contributions to 
the normal incidence given by state 1 (green, light gray) and state 2 (red, dark gray), for A = 4.0 and the same width of the 
barrier. (Lower Right) Contributions to the normal incidence given by state 1 (green, light gray) and state 2 (red, dark gray), 
for A = 4.0 and k = 5, for a varying width of the barrier. 


ii) k < A The system will change only through the 
transmitted state 1 within region II and reads 
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whose solution gives the TP plotted in Fig. 10. Here 
again we notice the already familiar behavior of selec¬ 
tion of states 2 only in the central region. We notice that 
the solutions behave just like a continuation of the fc > A 
problem and we have no energy resonances for the normal 
transmission problem. For small energy we have a trans¬ 
mission probability with poorly defined peaks, although 
if we look at the separate contributions of the two states, 
we see that they are such that the maxima and minima 
are opposite, except for normal transmission. Changing 
the size of the barrier, once again, induces resonances. 


VII. CONDUCTANCE AND THE LANDAUER 
FORMALISM 

We will try to use now the results obtained above to 
access a measurable quantity, namely, the conductance of 
the material due to the presence of the barriers. We un¬ 
derstand that our solution here might be an oversimplifi¬ 
cation, since we are not dealing with disorder, finite size 
effects, contacts, etc. Nevertheless, we believe that this 
might give us some insight on the usefulness of results 
like ours. We shall use a reasoning based on the Lan- 
dauer formalism. The total transmitted current through 
our sample should be given by 

j {E)Jl{E) 

(94) 

where 'D\ {E) is the one dimensional density of states per 
unit length for channel i, Jl {E) is the total transmitted 
current through our barrier through channel i and the 
sum is over all the channels, fii and ^2 are the respec¬ 
tive chemical potentials of two particle reservoirs thor¬ 
ough which the electrons will flow and /o is the Fermi 
distribution with null chemical potential. 

This leads to 


dfo {E) 
dE 


I _ I 
(Mi - M2) e A(() 


the dispersion relation is to be given by (putting back 
the fi’s) 


= p^=> E"^ =pI + {in/Lf 
pI = E^- {zir/Lf 
pi = E^-el 


where ei takes the value for py, pi the value for p^ and i 
are integer numbers. Notice that pi would give the energy 
for a Dirac particle moving in graphene in one dimension 
at the channel i. The total current we are dealing with 
here, Jl (E), is normal to the barrier, the possibilities of 
different incident angles is absorbed in the channels as 
we make clear below. 

Therefore, with this prescription, the transmission 
probability for the Vy interaction is given by 
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[ E "^ cos^ 0 — A^) 
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(96) 


where E (pi) = pi,which makes explicit our argument 
that the transmissions are ID-like, normal to the barrier, 
and separated in different channels. 

Now let’s analyze the expression for the current. We 
can write, 


J2eI{e)jI{e) = ^j;(F;)l^<j(F;-F;(p,)) 

i i Pi 


= \ Y.HE-E{p,))JI{E), (97) 

Pi,i 
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e I dEj2'E\{E)Jl{E) 


dfo (E) 
dE 


(95) 


where A(j) is the potential difference and G is the con¬ 
ductance between the reservoirs. 

The channels i are given by the discrete values that ky 
would assume in a finite sized sample, in such a way that 


and using the results of sec. IV for the transmission direc¬ 
tion only, Jl {E) = eT* {E) \vx {pi)\ = eT {E (p*)) \vx (pi)l 
( the transmission is assumed to be along the x direction 
and the modulus is present to guarantee that only the 
positive X direction is to be considered), one has 

J2e\{e)ji{e) = 
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Figure 10: (color online) Behavior of the TP for the Vx fc < A barrier problem . (Upper left) We have d — 3.0 and A = 6.0 
fixed and vary the energy k of the particles as 3.0 (green, full line light gray), 4.0 (red, mid gray) and 5.0 (blue, dark gray). In 
the middle, we separate the contributions to the total TP from the states l(full lines) and 2 (dashed lines) for the same values 
of d, A, k. (Lower left) The figure shows the contributions to the normal incidence given by state 1 (green, light gray) and 
state 2 (red, dark gray), for A = 4.0. Notice the continuation of this last graphic in comparison to the lower left graphic in 
figure 9. (Lower right) Contributions to the normal incidence given by state 1 (green, light gray) and state 2 (red, dark gray), 
for A = 4.0 and fc = 3 , for a varying width of the barrier. 


= - E{pi))T{E{pi))\v^{pi)\ 
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Pi 

= e^TiE)J26{E-E^{p^))\vAP^)\, (98) 

Pi 

where L is the length of the sample. Now we can stop 
looking at the values of Py as different modes and sum 
over all the values of the vector k ( in two dimensions) 
obtaining 


or, putting back the Fermi velocity(remember that in our 
units were such that vp = i) 

2W 

M{E) = —\E\. (101) 

Trnvp 

We can easily determine, now, the value of the conduc¬ 
tance due to our barrier. For zero temperature (in which 
case the integral is trivial) we get, 

2W 

G=—T{Ef)^—Ef ( 102 ) 
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(99) 


The function M (E) defines our number of transverse 
modes. It can be calculated by the usual method, trans¬ 
forming the sum in an integral and remembering to count 
the spin degeneracy. Valley degeneracy should not be in¬ 
cluded because the electrons are chosen in a very well 
defined valley, otherwise we should add another factor of 
2 multiplying the final expression. It gives, for a sample 
of width W 

2W 

M{E) = —\E\. (100) 


(notice that the units are correct). 

Hence, we show that the conductance is expected to be 
proportional to the TP in the normal direction with the 
energy set at Ep- As we have shown and analyzed in the 
last section the behavior of this TP we will neither bother 
writing the explicit expressions nor going through the 
analysis of the graphic representation once again. The 
only comment we believe it deserves is to notice that for 
the Vx case we expect to be able to choose the contri¬ 
butions to the conductance for each different cone, as a 
function of the width of the barrier. We also expect oscil¬ 
latory behavior of the conductance as we vary the energy 
in the case Vy. 


VIII. CONCLUSIONS 

We analyzed the problem of electronic transmis¬ 
sion in graphene through interfaces between regions in 
which quasi-particles belonging to different valleys (Dirac 
cones) interact or not. The relevant Hamiltonians we 
have employed are seen to be able to either create a gap 
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in the quasi-particle spectrum or shift the Dirac cones 
with respect to each other. In the latter case, hole states 
become available to positive energy incident particles and 
the Klein paradox arises. 

The behavior of the TP indicate that, for both Hamil¬ 
tonians, in barrier and step like problems, there is focal- 
ization of the incident beam. In barrier problems, due 
to the behavior of the wavefunctions in each different re¬ 
gion, we see that it should be possible to have situations 
in which we can enter the system with particles with mo¬ 
menta around one of the Dirac cones and come out with 
a superposition of electronic states about both cones. We 
could also act in the reverse way and filter states from 
one specific cone out of a general superposition of elec¬ 
tronic states involving different valleys. If systems in 
which the electron-electron interactions we proposed can 
be isolated, this physical phenomenon could be useful for 
the development of " valleytronics", without dealing with 
edge modes [201 HI] • 

Another possibility is to use this kind of transmission 
to entangle a pair of originally separable electronic states 
belonging to different cones. This would be very useful 
if one wishes to employ a graphene sheet in the develop¬ 
ment of quantum processors. 

In order to observe these effects, one should measure 
the contact resistance of small graphene samples placed 


on appropriate substrates which would induce the de¬ 
sired electron-electron interaction. Some considerations 
towards this were given in the last section, where we 
showed that the behavior of the conductance, at very 
low temperatures, is expected to be proportional to the 
TP normal to the given barrier. 

Finally, a word of caution about the Hamiltonians we 
have used to induce the coupling between different val¬ 
leys. Although we have appealed to general arguments to 
propose the phenomenological forms we have employed, 
we do not yet know of any microscopic mechanism to de¬ 
duce them. Nevertheless, we believe that they can indeed 
be obtained from a more microscopic approach and shall 
be investigating this possibility in the near future. 
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